Enhanced solid-state multi-spin metrology using dynamical decoupling 
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We use multi-pulse dynamical decoupling to increase the coherence lifetime (T2) of large numbers 
of nitrogen-vacancy (NV) electronic spins in room temperature diamond, thus enabling scalable 
applications of multi-spin quantum information processing and metrology. We realize an order- 
of-magnitude extension of the NV multi-spin T2 for diamond samples with widely differing spin 
environments. For samples with nitrogen impurity concentration < 1 ppm, we find T2 > 2 ms, 
comparable to the longest coherence time reported for single NV centers, and demonstrate a ten- 
fold enhancement in NV multi-spin sensing of AC magnetic fields. 



The negatively-charged nitrogen-vacancy (NV) color 
center in diamond possesses many useful properties — 
long electronic spin coherence times at room tempera- 
ture, an optical mechanism for initializing and detect- 
ing the spin state, and the ability to coherently ma- 
nipulate the spin using electron spin resonance (ESR) 
techniques — which make it a leading solid-state system 
for scalable applications in quantum information and 
metrology, such as sensitive detection of electric and 
magnetic fields with high spatial resolution jjHi] or in 
bulk [tj |T3| • Recently, dynamical decoupling techniques 
have been used to reduce the effective interaction of sin- 
gle NV spins with other spin impurities in the envi- 
ronment, enabling significant improvements in the NV 
single-spin coherence lifetime (T2) 11-13] and AC mag- 
netic field sensitivity 13, 3]. In this Letter, we demon- 
strate the successful application of dynamical decoupling 
to large numbers of NV spins (> 10 3 ) in room tem- 
perature diamond [3, [TR [l6j | . We employ multi-pulse 
Carr-Purcell-Meiboom-Gill (CPMG) and XY control se- 
to improve both the NV multi-spin T 2 
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and AC magnetic field sensitivity by an order of mag- 
nitude. We find similar relative improvements for di- 
amond samples with widely differing NV densities and 
spin impurity concentrations. For some samples, the NV 
multi-spin T2 is increased to > 2 ms, where it begins to 
be limited by NV spin-lattice relaxation (Ti ss 6 ms). 
These demonstrations of the utility of dynamical decou- 
pling for solid-state multi-spin systems pave the way for 
scalable applications of quantum information processing 
and metrology in a wide-range of architectures, includ- 
ing multiple NV centers in bulk diamond [3, [HI, 2D 
(thin-layer) arrays 0, [ljl QjJ |2(J, and diamond nano- 
structur es 12 ll |22| ; as well as p hosp horous donors in sili- 
con [23I. l24l] and quantum dots (ia j. 

The NV center is composed of a substitutional nitro- 
gen (N) impurity and a vacancy (V) on adjacent lattice 
sites in the diamond crystal (Fig. 1(a)). The electronic 
structure of the negatively charged state of the NV center 



FIG. 1: (a) Nitrogen vacancy (NV) color center in a diamond 
lattice. NV electronic spin decoherence is dominated by 13 C 
nuclear spin and N electronic spin impurities, (b) Energy 
level structure of negatively charged NV center, (c) Hahn 
Echo and n-pulse CPMG control sequences. Timing of AC 
magnetic field to be measured is shown in green. 



has a spin-triplet ground state, where the m s — ±1 levels 
are shifted from the m s = level by ~ 2.87 GHz due to 
the dipolar spin-spin interaction (Fig. 1(b)). Applica- 
tion of an external static magnetic field along the N-V 
axis Zeeman shifts the m s — ±1 levels and allows one 
to treat the m s = 0, m s = +1 spin manifold (for exam- 
ple) as an effective two-level system. The NV spin state 
can be initialized in the m s — state with above-band 
laser excitation, manipulated with resonant microwave 
(MW) pulses, and read out optically by measuring the 
spin-state-dependent fluorescence intensity in the phonon 
sideband. 

The NV spin environment (i.e., spin bath) is domi- 
nated by 13 C and N impurities, randomly distributed in 
the diamond crystal. These spin impurities create time- 
varying local magnetic fields at each NV spin, which can 
be approximated as an average local magnetic field that 
fluctuates on a timescale set by the mean interaction 
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between spins in the bath, inducing rapid dephasing of 
freely precessing NV spins on a timescale T 2 * ~ 1 /is for 
typical spin impurity concentrations. By applying a sin- 
gle resonant microwave it pulse to refocus the dephasing, 
the Halm Echo sequence (Fig. 1(c)) decouples NV spins 
from bath field fluctuations that are slow compared to 
the free precession time P, Q ■ Application of additional 
control pulses, as in n-pulse CPMG (CPMG-n) and XY 
sequences, have recently been shown to decouple single 
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NV spins from higher frequency bath fluctuations 

In the present study, we used a wide-held fluorescence 
microscope to perform dynamical decoupling and magne- 
tometry measurements on large numbers of NV centers 
in three diamond samples with different NV densities and 
spin impurity environments. A switched 3- Watt 532-nm 
laser provided optical excitation of NV centers within a 
10 /im-diameter cross-section of each sample. NV spin- 
state-dependent fluorescence was collected by a micro- 
scope objective and imaged onto a CCD array. Resonant 
microwave control pulses for coherent manipulation of 
the NV spin states were applied using a loop antenna de- 
signed to generate a homogeneous Bi field over the sam- 
ple detection volume. Applying a static field (B ~ 70 
Gauss) along one of the four diamond crystallographic 
axes selected approximately one quarter of the NV cen- 
ters to be resonant with the microwave pulses. Each 
diamond sample consisted of an NV-rich layer grown by 
chemical vapor deposition on a non-fluorescing diamond 
substrate, such that all collected fluorescence could be 
attributed to the NV-rich layer. 

Sample A (Apollo) had a 16 /im-thick NV-rich layer 
with NV concentration ~ 60 ppb (measured by NV flu- 
orescence intensity), N concentration ~ 100 ppm (mea- 
sured by secondary ion mass spectroscopy), and 1.1% 
natural abundance 13 C concentration. The high N con- 
centration dominated NV decoherence in this sample, 
limiting the measured Hahn Echo multi-spin coherence 
time to T 2 ~ 2 fis. We applied CPMG-n sequences and 
determined the NV multi-spin coherence time as a func- 
tion of the number of pulses, T 2 , from the 1/e decay 
of the spins' coherent evolution as a function of the to- 
tal CPMG-n evolution period (Fig. 1(c)). Representa- 
tive measurements of NV multi-spin coherence decay are 
shown in Fig. 2(a), with extended by a factor > 10 
for n = 128 (Fig. 2(d)). Furthermore, we found that T 2 (rl) 
exhibited a power-law dependence on n: T^ n ' cx n s : with 
s = 0.65 ± 0.02 for Sample A, which is consistent with 
the value s « 0.67 found recently for single NV centers in 



similarly nitrogen-rich diamond samples 11] . These re- 
sults demonstrate that inhomogeneities in the spin bath 
do not limit the effectiveness of dynamical decoupling 
for extending solid-state multi-spin coherence times by 
at least an order of magnitude. 

The NV-rich layers of Samples B (Apollo) and C (Ele- 
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FIG. 2: Measurements of NV multi-spin coherent evolution 
using n-pulse CPMG control sequences for three diamond 
samples of differing NV densities and spin impurity environ- 
ments: (a) NV ~ 60 ppb, N ~ 100 ppm, and 1.1% 13 C; (b) 
NV ~ 0.2 ppb, N ~ 0.1 ppm, and 1.1% 13 C; and (c) NV ~ 0.6 
ppb, N ~ 1 ppm, and 0.01% 13 C, where the solid lines denote 
fits to the decoherence envelope of the form exp [— (t/t 2 ) p ]- 
Note that the time axis of (b) is plotted with a linear scale due 
to the periodic collapses and revivals in the NV spin coher- 
ence of Sample B associated with 13 C Larmor precession [26rj . 
These collapses and revivals occur in samples where the 13 C 
abundance is high enough to contribute significantly to NV 
decoherence. (d) Scaling of measured NV multi-spin coher- 
ence times with the number n of CPMG pulses: T2 oc (n) s . 



ment Six) had much lower NV and N concentrations than 
in Sample A: NV ~ 0.2 ppb, N ~ 0.1 ppm for Sample 
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B and NV ~ 0.6 ppb, N ~ 1 ppm for Sample C. How- 
ever, while Sample B contained 1.1% natural abundance 
13 C, Sample C was isotopically engineered to reduce the 
13 C concentration to 0.01%. Despite these differences in 
spin impurity environments, we measured similar Hahn 
Echo NV multi-spin coherence times (T2 ~ 300 fis) for 
the two samples. By applying n-pulse CPMG dynamical 
decoupling sequences, we extended the multi-spin T^ 1 
to m 2 ms for both samples (Figs. 2(b) and 2(c)), which 
is comparable to the longest coherence time reported for 



dynamical decoupling applied to single NV centers [13 1 
and beginning to be limited by NV spin-lattice relaxation 
(Ti ~ 6 ms). 

For Samples B and C we also observed a power-law 



dependence for T, 



(n) 



(Fig. 2(d)), with lower scal- 



ing powers than for the N-rich Sample A: s = 0.42 ± 0.02 
for Sample B and s = 0.35 ± 0.01 for Sample C. These 
sample-dependent scaling powers suggest that multi- 
pulse dynamical decoupling control sequences can serve 
as spectral decomposition probes of spin bath dynam- 
ics [271 ] . 

We next applied dynamical decoupling to improve the 
sensitivity of NV multi-spin magnetometry. In a stan- 
dard AC magnetometry measurement utilizing a Hahn 
Echo sequence, an oscillating magnetic field, b(t) — 
b ac sin [(27t/t qc ) t + </>], induces a net phase accumulation 
of the NV spin coherence, which is maximized when the 
full time of the Hahn Echo sequence is equivalent to the 
period of the AC magnetic field (r ac ) and the phase offset 
<j) is such that the control pulses coincide with nodes in 
the magnetic field (Fig. 1(c)). Under these conditions, 
the field amplitude b ac can be extracted from the mea- 
surement of accumulated NV spin phase with optimum 
sensitivity, given approximately by [l|: 
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Here C is a parameter that encompasses the measure- 
ment contrast, optical collection efficiency, and number 
of NV spins contributing to the measurement. The con- 
trast is modified by NV decoherence over the course of 
the measurement, described phenomenologically by an 
exponential factor with power p. The value of p is found 
to be sample dependent, in the range of 1 to 2.5, and is 
related to the dynamics of the spin environment and to 
ensemble inhomogeneous broadening 27]. 

In an AC magnetometry measurement utilizing n-pulsc 
dynamical decoupling, the sensitivity is given approxi- 
mately by Eq. [T] with two modifications: (1) the mea- 



surement time is increased by r ac 



and (2) the 



NV multi-spin coherence time is extended by T 2 —> Tin 1 * 
The resulting sensitivity is given by: 



V(n) 



71 T) 



1 



: exp 



2T 2 



(2) 



Because the measurement time increases linearly with 
the number of control pulses n, whereas the coherence 
time increases sub-linearly, there is an optimum num- 
ber of pulses n opt that yields the most sensitive measure- 
ment of an AC magnetic field of period T ac given a set of 
sample-determined parameters: 
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For a given sample, all the parameters except r ac are 
fixed and we can simplify Eq. [3] to n opt oc (l/r ac ) (1 - s ~> . 
From this relationship, we see that at higher AC frequen- 
cies, more pulses are needed to reach the optimum sen- 
sitivity, which can be understood intuitively by realizing 
that the high frequency regime corresponds to short time 
intervals between control pulses during which time there 
is very little contrast lost due to decoherence (r ac <C T 2 ). 
More pulses increase the sensitivity by allowing for a 
longer measurement time and subsequently more phase 
accumulation per measurement. This intuition also illus- 
trates why multi-pulse sequences are more effective at en- 
hancing magnetometry sensitivity in the high-frequency 
regime, where the Hahn Echo scheme provides relatively 
poor magnetic field sensitivity (Fig. 3(c)). Note that 
extension of the NV multi-spin coherence time via multi- 
pulse dynamical decoupling (and thus enhancement of 
magnetic field sensitivity) is eventually limited by NV 
spin-lattice relaxation (Ti), beyond which increasing the 
number of control pulses is ineffective. 

We used Sample C to perform NV multi-spin magne- 
tometry measurements comparing Hahn Echo and multi- 
pulse dynamical decoupling schemes. The diamond de- 
tection volume was approximately 30 /1m 3 , correspond- 
ing to ~ 10 3 sensing NV centers aligned along the static 
magnetic field. We employed rt-pulse XY sequences, in 
which control pulses are applied with the same timing 
as in CPMG-n sequences but with alternating 90° spin 
rotation axes to provide more isotropic compensation for 
pulse errors [11]. Fig. 3(a) illustrates how we extracted 
the magnetic field sensitivity from the measured vari- 
ation of NV fluorescence as a function of applied AC 
magnetic field amplitude [2, 0] for different measurement 
pulse sequences. The AC magnetic field sensitivity im- 
proved with the number of pulses, in good agreement 
with predicted values up to n ~ 150 pulses, at which 
point pulse infidelities began to degrade the measured 
sensitivity (Fig. 3(b)). 

Over a wide range of AC magnetic field frequencies 
our NV multi-spin measurements confirmed that multi- 
pulse dynamical decoupling outperformed the Hahn Echo 
scheme, in agreement with theoretical expectations (Fig. 
3(c)). The enhancement in magnetic field sensitivity pro- 
vided by multi-pulse dynamical decoupling was especially 
pronounced at frequencies higher than the Hahn Echo 
I/T2 coherence. For example, at a frequency of 220 kHz, 
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FIG. 3: Measured AC magnetic field sensitivity for a 30 (im 3 
sensing volume of Sample C (~ 10 3 sensing NV spins), (a) Ex- 
amples of measured normalized fluorescence signals as func- 
tions of AC field magnitude b ac using a Hahn Echo sequence 
(green) (wider amplitude range shown in the inset) and multi- 
pulse XY sequences with 28 (red) and 54 (blue) control pulses, 
illustrating how the uncertainty in the measured signal (8S) 
limits the uncertainty in the extracted magnetic field magni- 
tude (SB). The sine behavior of the signal with respect to 
b ac is achieved by shifting the phase of the last microwave 
pulse by 90° from what is shown in Fig. 1(c). (b) Sensi- 
tivity as a function of number of pulses for an 84.5 kHz AC 
field (circles). The solid line is calculated using Eq. [2] with 
experimental parameters C = 0.1, T2 = 250 /-is, p = 1, and 
s — 0.37. (c) Comparison of calculated (lines) and measured 
(points) sensitivity at several AC magnetic field frequencies. 



we measured a factor of 10 improvement in magnetic field 
sensitivity: from 67.7 ± S.SnT/v'Hz using a Hahn Echo 
sequence, to 6.8 ± OAnT/i/Hz using a 240-pulse XY se- 
quence. 

In summary, we experimentally demonstrated that 
multi-pulse dynamical decoupling sequences extend by an 



order of magnitude the coherence lifetime of large num- 
bers (> 10 3 ) of NV electronic spins in room temperature 
diamond, for samples with widely-differing NV densities 
and spin environments. We realized an extension of the 
NV multi-spin coherence time to > 2 ms, which is com- 
parable to the best results from application of dynamical 
decoupling to single NV centers. We also showed that 
multi-pulse dynamical decoupling improves NV multi- 
spin AC magnetic field sensitivity relative to the Hahn 
Echo scheme, with a ten- fold enhancement for higher 
frequency fields. Further improvements in NV multi- 
spin coherence time and magnetic field sensitivity are ex- 
pected from the integration of dynamical decoupling with 
diamond samples engineered to have optimized spin en- 
vironments, as well as quantum-assisted techniques [29j 

. The 



and increased optical collection efficiency |21|, [ 
present results aid the development of scalable applica- 
tions of quantum information and metrology in a variety 
of solid-state multi-spin systems, including NV centers in 
bulk diamond, 2D arrays, and nanostructures; phospho- 
rous donors in silicon; and quantum dots. 
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